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Aims. We present a technique to determine the orbital and physical parameters of eclipsing eccentric Wolf-Rayet + 0-star binaries, 
where one eclipse is produced by the absorption of the O-star light by the stellar wind of the W-R star. 

Methods. Our method is based on the use of the empirical moments of the light curve that are integral transforms evaluated from 
■ the observed light curves. The optical depth along the line of sight and the limb darkening of the W-R star are modelled by simple 

' mathematical functions, and we derive analytical expressions for the moments of the light curve as a function of the orbital parameters 

and the key parameters of the transparency and limb-darkening functions. These analytical expressions are then inverted in order to 
(~ i , derive the values of the orbital inclination, the stellar radii, the fractional luminosities, and the parameters of the wind transparency 

^ |. and limb-darkening laws. 

I ' Results. The method is applied to the SMC W-R eclipsing binary HD 5980, a remarkable object that underwent an LBV-like event in 

O ' August 1994. The analysis refers to the pre-outburst observational data. A synthetic light curve based on the elements derived for the 

system allows a quality assessment of the results obtained. 
C/3 . 
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> 1. Introduction 

Photometric monitoring of Wolf-Rayet (W-R) binaries revealed 
\^ , that many of them display a shallow eclipse when the W-R star 
' passes in front of its O-type companion (e.g. Lamontagne et al. 
, — T [ 11996 ). These so-called atmospheric eclipses arise when part of 
. the light of the O-type companion is absorbed by the wind of the 
W-R star In a few cases, the light curve displays an eclipse at 
both conjunctions and the analysis of this phenomenon can pro- 
• • vide important information about the physical parameters of W- 
^ R stars and their winds. In this context, the most famous example 
k> J is V444 Cyg (WN5 + 06 V), which has been extensively inves- 
j_] ■ tigated by the Moscow group (e.g. Antokhin & Cherepashchuk 
d ' 120011 and references therein). Cherepashchuk and coworkers 
' (e.g. Cherepashchuk 119751 Antokhin & Cherepashchuk 120011 
and references therein) developed a sophisticated method to han- 
dle the ill-posed problem of light curve inversion for V444 Cyg. 
Based on the minimum a priori assumptions about the trans- 
parency function, this method not only yields the radii of both 
components, but also provides information about the structure of 
the WN5's stellar wind. However, because of a number of funda- 
mental hypotheses that are not necessarily valid for all eclipsing 
W-R binaries, this method cannot be readily applied to all eclips- 
ing W-R + O systems. For eccentric systems in particular (e.g. 
WR22, Gosset et al. 119911 ) some assumptions (such as spherical 
symmetry of the problem) break down and a different technique 
must be used. 

We initiated our study of the observed light changes of 
eclipsing Wolf-Rayet binary systems when two of us (J.B. and 
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CP.) tried to confirm the 25.56 day period found by Hoffmann 
et al. (fT978l) for the SMC star HD 5980, the first extragalactic 
Wolf-Rayet binary then known to display eclipses. This exercise 
led to the discovery of the correct orbital period of HD 5980, 
P = 19.266 + 0.003 days (Breysacher & Perrier fTMOl l, the light 
curve revealing, in addition, a rather eccentric orbit of e = 0.47 
assuming / = 80°. However, because of the uncertainties in the 
depth of both minima, caused by an insufficient number of ob- 
servations, no detailed quantitative analysis of this preliminary 
light curve could be attempted. 

Its relatively long period and large eccentricity ensure that 
HD 5980 is an interesting object in which to study the structure 
of a W-R envelope, and a photometric monitoring of this sys- 
tem was initiated to define the shape of its light curve in a more 
accurate way. More than 700 observations were collected. After 
realizing that none of the existing 'classical tools' was suited to 
our purpose - the decoding of the light changes of a partially- 
eclipsing system characterized by an eccentric orbit and contain- 
ing one component with an extended atmosphere - we started to 
develop another approach to the solution of light curves. 

The technique of light curve analysis applied to V444 Cyg 
by Smith & Theokas (1980), which is based on Kopal's fun- 
damental work (cf. Kopal , 1 9751 1 1 9791 1, appeared as an attractive 
approach to the solution of our problem. This method is based on 
the interpretation of the observed light changes in the frequency- 
domain, i.e., not the light curve as a function of time, but its 
Fourier-like integral transform. 

We now describe in detail the method we developed for the 
study of Wolf-Rayet eclipsing binaries. A preliminary applica- 
tion of this technique to the light curve of HD 5980 prior to 
its 1994 LBV-like eruption (see e.g., Bateson & Jones 119941 
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Barba et al. 119951 Heydari-Malayeri et al. 1 19971 1 was presented 
by Breysacher & Perrier (1991 hereafter BP91). We reanalyse 
the pre-outburst observational data using an improved version 
of the software tool. Revised values for the physical parameters 
of HD 5980 are derived. A synthetic light curve based on the 
elements thus obtained allows a quality assessment of the new 
results. 

A more detailed study of HD 5980, including the analysis of 
the light curve obtained after the eruption (Sterken & Breysacher 
119971 ). will be presented in a forthcoming paper. 

2. Analysis of the light changes in the 
frequency-domain 

In this section, we present the fundamental equations of our 
method, we then introduce the mathematical functions to model 
the transparency and limb-darkening functions and consider the 
specific problem of eccentric orbits. 

2.1. The basic equations 

We refer to the fundamental work of the Manchester group (cf. 
Kopal [T975l 119791 Smith [T976|, and first consider an eclipsing 
system consisting of two spherical stars revolving around the 
common centre of gravity in circular orbits, and appearing in 
projection on the sky as uniformly bright discs. The system is 
seen at an inclination angle /. When star 1 of fractional lumi- 
nosity Li and radius ri is partly eclipsed by star 2 of fractional 
luminosity L2 and radius r2 (Fig.[T]), the apparent brightness / of 
the system (maximum light between minima taken as unity) is 
given by 

l(rur2,S,J)^ 1 - J J J(r)dcr, (1) 

where 6 is the apparent separation of the centres of the two discs, 
J represents the distribution of brightness over the apparent disc 
of the star undergoing eclipse, and dcr stands for the surface el- 
ement. The distances ri, r2, and 6 are expressed in units of the 
orbital separation. The integral in Eq. ([TJ provides the apparent 
Toss of light' displayed by the system when an area A(ri, r2, 5) 
of star 1 is eclipsed (see Fig.[TJ. The assumption that star 1 is 
uniformly bright yields 

J(r) = ^. (2) 
Combining Eqs. (HI) and ©, we obtain 

1-Z(ri,r2,5,7)= f f do- = all, (3) 
7rr| J Ja 

where a is the ratio of the mutual area of eclipse to the area of the 
disc of the eclipsed star, and is a function of r^, r2, and 6 (Kopal 
11975 ). A generalisation of these concepts to the case of spherical 
stars with arbitrary limb-darkening laws J(r) was presented by 
Smith ( IT976b . 

For an orbital period P and an epoch of conjunction to, we 
define the phase angle at a time t to be 

2n 

e^—it-to). (4) 

As proposed by Kopal ( 119751 119791 ), we focus our attention on 
the area subtended by the light curve in the (I, sin^'" 0) plane. 



Star 2 




S 



Fig. 1. Geometry of the eclipse of star 1 by star 2. The integral of 
Eq. ([T) is evaluated over the hatched area A(ri,r2, 6). The points 
inside this area can be specified by the coordinates (r, s) corre- 
sponding to the two intersections of the circle with radius r and 
centre Oi and the circle with radius s centred on O2. 



where m is a positive integer number (m — 1, 2, 3, ...), as shown 
in Fig.|2] The areas A2m between the lines / = 1, sin^"' 6-0, and 
the true light curve are then given by the integrals 

A2„, = {l-l)d(sm^"'0), (5) 
Jo 

which are hereafter referred to as the moments of the eclipse, of 
index m, where 0^ denotes the phase angle of the first contact 
{6(6 fc) = ri -H r2) of the eclipse. 




Fig. 2. Light curve of an eclipse in the {I, sin " 6) plane. 6fc corre- 
sponds to the phase angle of first contact, whilst di represents the 
phase angle corresponding to the beginning of the total eclipse. 
The shaded area illustrates the moment A2m and A - 1-1(6 - 0). 

Combining Eqs. (□ and Kopal (1 19751 19791 ) and Smith 

(119761 ) demonstrated that, based on certain assumptions, it is 
possible to 
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• derive analytical expressions of the moments of the eclipse in 
terms of the physical parameters (/, r\,r2,L\, L2) of a binary 
system with a circular orbit consisting of uniformly bright 
spherical stars, 

• invert these relationships to determine the parameters of the 
system in terms of the moments A2,n that can be empirically 
obtained from the data, 

• extend this treatment to the case of partial eclipses of stars 
with an arbitrary (yet analytical) limb-darkening law. 

Going one step further. Smith & Theokas ( 1980 ) generalised 
the above concepts to derive a convenient mathematical solu- 
tion to the problem of an atmospheric eclipse, i.e., an eclipse 
of a limb-darkened star by a star surrounded by an extended at- 
mosphere. In such an eclipse, at each position inside the area A 
specified by the coordinates r and s (see Fig.[T]), a fraction of the 
light emitted by star 1 is absorbed by the atmosphere of star 2. 
To account for these transparency effects, a transparency func- 
tion F{s) is introduced into Eq. ([TJ, so that the total amount of 
light seen by the observer becomes 



lirur2,6,J,F) 



-/I 



J{r)F(s) do-. 



(6) 



Considering that in the case of an atmospheric eclipse, it might 
be interesting to give more weight to the data close to mid- 
minimum. Smith & Theokas (1 19801 1 also introduced an alterna- 
tive set of moments B2m defined by 



B2n 



re, 
Jo 



(1 



dicos^"'0). 



(7) 



The use of kernel d(cos^"' 9) in Eq. O places more emphasis on 
the data points close to mid-eclipse, which leads to smaller er- 
rors than in the case of the A2,,, moments defined by means of 
the t/(sin^"' 6) kernel (Theokas & Smith 1983, private communi- 
cation). We define e to be the mean error in an individual data 
point; the relative error in the light curve s/(l-l{0)) increases for 
data points near Of^ . While these points are given more weight by 
the kernel of the A2m moments, the converse situation holds for 
the B2m moments, where the kernel reaches its peak for a given 
m in a zone where 1 - 1(0) is closer to its maximum value. 

By definition, the A2m and Z?2m moments are related to each 
other by means of 



B21 



=z 



(m - p)\ p\ 



-{-\)P^^A2p. 



(8) 



We now concentrate on the 82,,, moments because they are 
equally well suited to the analysis of the primary and sec- 
ondary minima with the transparency and limb-darkening func- 
tions adopted in the present study (cf. Section 2.2). Since there 
are a number of typos in the paper of Smith & Theokas (1 19801 1, 
we provide below the mathematical details of the method. 

The infinitesimal element of area dcr of Eq. (|6]l is expressed 

as 



dcr - 



1 52 



2drd^ 



■ {nP'air, s, 5)) drds. 



(9) 



The following general expression for the B2m moments was de- 
rived by Smith & Theokas a980i l 



B2n 



rr 

Jo Jo 



J(r) F{s) dsdr 



drds 



2,n^d(nr^-a(r,s,6))^^ 
do 



(10) 



where 1 - 1(0 - 0) is defined as A (see Fig.|2]i. 

The analytical expressions obtained for the B2m's for m = 
1, 2, 3, 4, and 5 are thus 



B2 ^ A-csc^i(P- IiRir2-tffi), 



(11) 



84 = A- esc"* i(P -lIiRirl + l2R\r\rl + hRir^ - 1^2), (12) 



Be = A-c^.e'iiP-^hRi^+'^hRxrA + 'ihRirl 
-hRir\rl - hR^rl - 3 hRir]^ - iffi). 



A - csc^ i (P - 4 hRirl + 6 hRy^rl + 6 11R24 
-4 hRir\rl - 4 /i/?3r^ - 12 hRir^rj + ^4^1 '"1 '"2 
+IiR4rl + 6 hR2rtrt + 6 hR^rj^ ' 'A4), 



(13) 



(14) 



and 



Bio = A - esc'" i(P- 5 IiRi4 + 10 hRir^iij 
+ lQhR2r2 - 10 hRir'^^rl - IQIiRjr^ 
-30l2R2r^^4 + 5 hRir'irl + 30 hR2r'l4 
+30 hR^r^r^ + 5 IiR4rl - hRir\rl 
-\OI4R2r\r\ - 20 hRir\r\ - \OhR4r\rl 
-hRsA" - ilfs). 



(15) 



The coefficients P, /,„, /?,„, and tpm are defined (cf. Smith & 
Theokas 1980) by the equations 



P(rur2,J,F) 
hn{r\,J) 



•min{r\ ,r2) 



J{r)F{r)2nrdr, 



Rm(r2,F) 
and 



r 

Jo 

/ 

Jo 



'^'^ \nr''")dr. 



„2m-2 



dr 



ns) d 



r|™ ds 



(s^'")ds. 



(16) 
(17) 

(18) 



lAm = tlf„(ri,r2,i,J,F) = 



\ +2m(ri + r2f 
\ 6 (ri + r2Y 



x(cos^ ;■ - (r2 - rif-)^^^ - Vcos^ ; - (r2 - ri)^ 
m\r2 - nl 
8(ri +r2)2 



(cos'*/-(r2 -ryf) 



+\r2 — ri\ arctan 



-y/cos^ ; - (r2 - ri)2 



K2 - '-I I 



(19) 



2.2. The transparency and limb-darkening functions 

In our method, the transparency of the W-R wind is described 
by an analytical function that depends on a limited number of 
parameters. Since the functional form of the transparency is 
adopted a priori, our choice will obviously have a direct influ- 
ence on the parameters derived for the system. Therefore, it is 
important to clearly specify the assumptions made in our ap- 
proach. To avoid confusion with the standard symbols employed 
by Smith and Theokas (1980), from now on, the various radii 
in our model will be denoted p, (i=l,2,3), where pi stands for 
the radius of the O-star that undergoes the eclipse. The use of 
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the subscripts e and a will indicate whether the components are 
seen in emission or absorption. 

The first simplifying hypothesis is that all composite parts 
of the system are supposed to be spherically symmetrical. This 
means that the method is not applicable to close binary systems 
in which the components depart strongly from a spherical form 
as a result of tidal distortion and where ellipticity and reflection 
effects are both present. 




1-y 



y 



..Fj (opaque W-R core) 



(semi — transparent W— R envelope) 
V =0 




Fig. 3. Schematical view of the transparency law across the disc 
of the W-R component as given by Eq. (ISTT i. The solid line shows 
the amount of light absorbed along the line of sight of impact 
parameter s. The individual contributions due to the opaque core 
and the semi-transparent extended atmosphere are illustrated. 

Because a principle objective of the proposed technique of 
light curve analysis is the determination of the structure of a W- 
R envelope, a composite model consisting of an opaque core 
and a surrounding extended atmosphere was adopted for the 
W-R component. As a consequence, while for the transparency 
function F{s) of the eclipsing W-R star, of radius ro. Smith & 
Theokas (,1980) simply adopted 



F(s) = ^,,,('"0'^) = Va 



l-v{- 



for s < ro. 



(20) 



where v is the coefficient of transparency, we chose as a first step 
a transparency law of the form 



F(i) = F,_,,>3a, 0) + Fv>2«, l^), 



(21) 



where the radius of the opaque core of the W-R star is p^a and 
that of the extended eclipsing envelope is p2a (> p3a, see Fig. [3]). 
This transparency law was used in the preliminary analysis of 
the light curve of HD 5980 by BP91 and is motivated by the fact 
that it corresponds to a physically more realistic model of the 
W-R star than that defined in Eq. ( |20| ). although the advantage of 
relative mathematical simplicity is still preserved. 

For the brightness distribution J{r) over the W-R disc (im- 
portant for the eclipse of the W-R component by the O-star), a 
law very similar to that of the transparency function is adopted 



J(r) = 7(0) [jl-,Me, 0) + JyMe, «2)] , 



where 7(0) is the central surface brightness, U2 is the coefficient 
of limb-darkening, and 7,^ is defined as 



-/>v('"0,m) = ye 



l-u\ — 



(23) 



When the W-R star is eclipsed, the radius of the core, assumed 
to be of uniform brightness, becomes and that of the limb- 
darkened envelope p2c- We note that the core and envelope radii 
of the W-R component seen in emission or absorption may differ 

We briefly address the physical meaning of this transparency 
law. First, in our model, the semi-transparent envelope of the 
W-R star has a finite extension given by the radii p2a and p2e- 
However, the stellar winds of W-R stars do not stop abruptly 
so close to the star, but instead extend to large distances (much 
larger than the orbital separation in most binaries) to the shock 
with the interstellar (or circumstellar) medium. How can we then 
interpret p2a and p2c? The radius p2a corresponds to the farthest 
position in the stellar wind where the residual optical depth along 
the line of sight produces a variation in the light curve that can be 
distinguished against the intrinsic photometric variability of the 
W-R star and the photometric errors. Similarly p2e is the outer 
radius of the W-R envelope that emits a measurable fraction of 
the light in the considered waveband. 

A clear difference between our approach and that of 
Antokhin & Cherepashchuk (2001) concerns the functional be- 
haviour of F(sy. whereas in our model, F{s) is a convex function 
over the entire range p e [0,p3a], Antokhin & Cherepashchuk 
(2001) use a convexo-concave function, where the concave part 
corresponds to the stellar wind. As a consequence, F{s) given by 
Eq. in} decreases at a slower rate over the wind than the trans- 
parency law inferred by Antokhin & Cherepashchuk. 

In contrast to Smith and Theokas (19801 1 who neglect the ef- 
fect, we take into account the limb-darkening of the OB-type 
component. Assuming that the formula employed by these au- 
thors to represent the brightness distribution across the W-R 
disc also applies to normal stars, we adopt the following limb- 
darkening law for the OB star 



Lo 



;rpj(l - ui + Mj/3) 



1 - Ml I — 

Pi 



(24) 



where Lq is the luminosity of the OB star, of radius pi, and ui is 
the coefficient of limb-darkening at the effective wavelength of 
the photometric filter considered. 

Using these laws of transparency and limb-darkening, we 
then derived the expressions for P, /,„, (see AppendixlAli. and 
tpm, and hence the final equations for the moments B2m, corre- 
sponding to the primary and secondary minima. 

2.3. The orbital eccentricity 

The treatment of elliptical orbits in the frequency-domain was 
also addressed by Kopal (I1979I I. The problem still concerns the 
determination of the elements of the eclipse from the moments 
- B2m in the present case - derived from the light curve, but 
accounting for the eccentricity e and the longitude of periastron 

0). 

In the definition of the B2m moments, the phase-angle 6 is 
no longer identical to the mean anomaly M but has rather to be 
replaced by a linear function of the true anomaly v 



(22) 



7T 

V + a> . 

2 



(25) 
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The (i(cos^™ 9) kernel in Eq. (|7J thus becomes 



c/(cos2'" 0) 



d [sin^'" (v + u))\ . 



(26) 



As a consequence, the empirical values of B2m cannot be derived 
from the observed data until a proper conversion of the phase 
angle into the true anomalies has been completed. This can be 
achieved either by a numerical inversion of Kepler's equation or 
the well-known asymptotic expansion of elliptical motion (e.g. 
Danion ri959. Kopal & Al-Naimiv [T978b 



V = M + (2e- ie^)sinM + (-e^- — e'*)sin2M 
^ 4 ^ U 24 ^ 

13 , 103 4 

H e sm3 M H e sm4M + ... 

2 96 



(27) 



Regardless of the technique used to compute the true 
anomaly, this conversion evidently requires an a priori knowl- 
edge of e as well as w. For a given value of the inclination /, 
these parameters can be derived by inversion of the equations 
(see e.g. Kopal & Al-Naimiv [19781 l 



AO 



and 



1 



1 1 + CSC / 



— [- cos^ w - 2 + C>(cot'^ /)]), 



e sm Lo — 



d2 - di 



1 - 



sin^ 



(28) 



(29) 



where AO, di, and dx are, respectively, the phase displacement of 
the minima and the durations of the primary and secondary stel- 
lar core eclipses. These latter quantities are determined directly 
from the observed light curve. Since the orbital inclination of an 
eclipsing binary system is likely to be rather large, the (9(cot^ /) 
term in the coefficient of the term of Eq. ( l28l l can be neglected. 

The empirical 'elliptical' moments of the light curve then 
provide the elements of the binary exactly as in the 'circular' 
case. One must ensure that the resulting values of the radii have 
been reduced to a constant unit of length. In our code, we there- 
fore report all distances in relation to the semi-major axis a of 
the relative orbit. 

Our composite model adopted for the W-R star has differ- 
ent radii {p2,3a and p2,3e) depending on whether the W-R com- 
ponent is seen as an eclipsing or an eclipsed disc. Because of 
this, although each individual half-eclipse provides an indepen- 
dent solution for the elements, the complete determination of the 
elements requires a combination of solutions obtained for both 
minima and because of the non-zero eccentricity, both the de- 
scending and ascending branches of each. 

3. Decoding of the light curve 

3.1. Empirical determination of tiie moments 

We consider the light curve of a W-R binary system of ec- 
centricity e and period P, derived for a given photometric 
bandpass. The data are assumed to consist of a list of entries 
that provide for each observation the orbital phase O, and 
the measured intensity Z,. To normalize the brightness scale, 
a mean intensity value is derived well outside the eclipses, 
during a phase-interval where the system is assumed to display 
(constant) maximum light. The /,■ value of each data point is 



then divided by this mean to normalize the light curve to unity. 
It has to be emphasized, however, that this does not necessarily 
imply that Li + L2 = 1 for the W-R binary. The luminosity of 
a third photometrically unresolved component along the line of 
sight may indeed contribute to the observed brightness as well, 
thereby leading to a brightness distribution such as Li -H L2 < 1 
in the final solution for the eclipse. In the case of an eccentric 
orbit, the eclipse-free mean Z, value is preferably taken around 
apastron to avoid as much as possible any luminosity increase 
that could occur around periastron as a result of enhanced 
interaction effects between the components. 

The determination of the moments B2m requires a smoothed 
light curve, which can be obtained, for instance, from a spline 
fit to the observed points with special attention to the minima. 
However, this task can become difficult if the descending or as- 
cending branch of either minimum is ill-defined because of an 
uneven sampling of the observations or intrinsic photometric 
variability in the W-R star (see e.g. the case of WR22, Gosset 
et al. 1199 1'). A clustering of the points, in particular, is a serious 
handicap for the method. A preliminary processing performed, 
by filtering the observational data, to help reduce the scatter al- 
lows us then to obtain a smooth light curve Z(0). 

The next step consists of determining the quantities di, d2, 
and AO (see above) from the smoothed light curve. For an as- 
sumed value of the orbital inclination /, the parameters e and a> 
are obtained by means of an inversion of Eqs. ( l28T l and ( |29] ). With 
these values of e and co, the orbital phases O, are converted into 
true anomalies. 

The moments 82,,,, which take into account the eccentricity 
effect, are obtained in practice by summation, using the follow- 
ing expression 

N-l 

B2m = 2(COS2'«(0,.)-COS2'«(0,.+ 1)) 



(30) 



where 0, are the predefined angles at which the smoothed light 
curve is sampled, is defined by the constant step A0 - 0,+ i - 
0, adopted, and the value of ©i corresponds to the first contact 
of the eclipse. The Z(0,) values refer to the normalized smoothed 
light curve, and by definition I - 1 for |0| > |0i|. 

Since the individual data points are affected by observa- 
tional errors, the integration of the empirical moments must itself 
be affected by errors. The uncertainty associated with the mo- 
ments B2m can be evaluated using the following equation (see 
Al-Naimiv, [T977l Smith & Theokas [T980l l 



AB 



2m 



X (cos^'" 61 - cos^™ 9„) , 



1/2 



(31) 



where n is the number of observed points over the considered 
eclipse, and Ij -l{0j) is the difference between the observed point 
of index j and the smoothed light curve at 6j. The angles 0\ and 
0„ refer, respectively, to the first and last observed data point over 
the relevant part of the light curve. 

3.2. Solution for the elements 

For each half-eclipse, there are five non-linear algebraic equa- 
tions to be solved simultaneously for the elements pi, p2fl or p2e. 
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p-ia or p3e, /, Ml or U2, ja oTje, V, Li ov Lo the meanings of which 
are summarized below for convenience: 

Li - luminosity of the OB-type star, 

L2 = luminosity of the W-R star, 

i = inclination angle of the orbit, 

pi - radius of the OB-type star, 

P2a,e = radius of the W-R envelope seen in absorption or 
emission, 

P3a,e - radius of the W-R opaque core seen in absorption or 
emission. 

Ml = limb-darkening coefficient of the OB-type star, 

M2 = limb-darkening coefficient of the W-R envelope, 

u = transparency coefficient of the W-R envelope (cf. Fig. [3]), 

and 

ya,e - contribution of the W-R envelope in absorption or 
emission (cf. Fig. [3]) 

This large number of variables can fortunately always be re- 
duced to a smaller number for both the primary and secondary 
minima, as explained hereafter The preliminary determination 
of the orbit inclination / already eliminates for instance one vari- 
able. 

At the primary minimum, when the OB star is in front, ac- 
cording to our composite model v = Qhy definition. Since we 
must also have that < M2 < 1, solutions can be searched for a 
set of discrete values of the parameter M2 in this interval, so that 
the remaining variables are p\, p2e, P3e, ye, and L2- 

At the secondary minimum, when the W-R star eclipses the 
OB component, mi is the limb-darkening coefficient of the OB 
star. The value of mi can be adopted following e.g. the tabulated 
values supplied by Klinglesmith & Sobieski d 19701 1. According 
to our model, we now have < u < 1, so that again, after se- 
lection of a sample of v values, we can proceed in solving the 
equations for the remaining parameters pi, p2a, p3a, ya, and Li 
only. 

The solution of the system of as many as five non-linear 
equations 

B2m(pi,P2a,Pia,i,ui,ya,v,Li,L2) = B2m(observed) 

or 

B2m(pi,P2e,pie,i, U2,ye,v, U, L2) = B2m(observed) 

is obtained by minimizing the compiled from the residuals of 
these equations 

N 2 

2 V |B2m(computed) - B2m(observed)r 

^ - L xbI ■ ^'^^ 

m=l 2m 

This minimization is achieved by means of Powell's technique 
(e.g.. Press et al. 1 19921 and references therein). 

4. Application to HD5980 

4.1. HD 5980: a peculiar system 

HD5980 s AB5 (Azzopardi & Breysacher'1979) is associated 
with NGC 346, the largest H n region + OB star cluster in the 
Small Magellanic Cloud. This remarkable W-R binary, which 
underwent an LBV-type event in August 1994, is presently rec- 
ognized as a key-object for improving our understanding of mas- 
sive star evolution. HD 5980 is a rather complex system because 
it consists of at least three stars: two stars form the eclipsing bi- 
nary with the 19.266 day period, whilst the third component, an 
O-star, which is detected by means of a set of absorption lines 



and by means of its third light (see also below), could be a mem- 
ber of a highly eccentric 96.5 day period binary (Schweickhardt 
|2000 Foellmi et al. 2008). Whether or not the third star is phys- 
ically bound to the eclipsing binary remains currently unclear 
Before the LBV eruption, both components of the eclipsing bi- 
nary akeady showed emission lines in their spectra and were 
thus classified as Wolf-Rayet stars (Niemela[1988jl. However, as 
shown by the analysis of the spectra taken during and after the 
LBV event, at least the star that underwent the eruption (here- 
after called star A) was not a classical, helium-burning, Wolf- 
Rayet object, but rather a WNha star, i.e., a rather massive star 
with substantial amounts of hydrogen present in its outer lay- 
ers (Foellmi et al. 2008 ). These WNha stars have wind proper- 
ties that are intermediate between those of extreme Of stars and 
classical WN stars. 

A summary of the light changes exhibited by HD 5980 was 
presented by Breysacher (I1997I I. The technique of light curve 
analysis described above is applied to HD 5980 prior to the out- 
burst. Given that star A, the component in front of its companion 
(hereafter called star B) during primary eclipse, was a WNha star 
and since there are no indications of wind effects in the primary 
eclipse, we assume that this component behaves as an OB-star in 
the light curve. 705 measurements with the Stroemgren v filter, 
described in BP91, are taken into consideration. The shape of 
the resulting light curve does not allow us to use the ill-defined 
ascending branch of the primary eclipse (star A in front) for the 
analysis. Therefore, only three half-minima will be considered. 
Compared to the preliminary analysis carried out by BP91, the 
software tool presently used has been upgraded, allowing us for 
instance to assess the quality of the solution by means of a syn- 
thetic light curve. 

4.2. Solutions of the light curve 

From the smoothed light curve, the durations of the primary and 
secondary stellar core eclipses as well as the separation between 
the core eclipses are measured first to be di - 0.062 + 0.005, 
d2 = 0.095 ± 0.005, and AO = 0.362. All durations are ex- 
pressed as phase intervals (i.e. fractions of the orbital cycle). The 
corresponding values of the eccentricity e and the longitude of 
periastron w are: e = 0.314 + 0.007 and w = 132.5°+ 1.5°. These 
values are in fairly good agreement with those derived for these 
parameters by Breysacher & Francois ( I2000I I (e - 0.30 ± 0.02, 
a> = 135° ± 10°) from a completely different approach based on 
the analysis of the width variation of the He ii A 4686 line using 
the analytical colliding- wind model devised by Liihrs ( I1997I I. 
From radial-velocity studies, Kaufer et al. (120021 and Niemela 
et al. ( [19971 ) also found that e = 0.297 + 0.036 and e = 0.28, 
respectively. The values of the moments of the light curve of 
HD 5980 are listed in Tableffl 

The inclination angle of the orbit can easily be derived with 
reasonable accuracy. For each of the three half-minima, a quick 
analysis is carried out for a number of plausible discrete val- 
ues of / (/ = 82°, 83°, ...89°), and the value of / finally adopted 
is the one for which closest agreement is obtained between the 
three solutions provided for the radius of star A and the radius 
of the opaque core of star B. These quantities are fundamental 
elements of the system and the combination of the solutions of 
the descending and ascending branches of both minima must in- 
deed provide, at the end of the detailed analysis, a unique value 
for pi and very similar - if not identical - values for p3a and p3e. 
The above conditions are fullfilled for 85° < / < 87°, therefore 
we adopt / = 86° ± 1°. 
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Table 1. Values of the moments of the eclipses of HD 5980 used 
throughout this paper. These values are derived from the light 
curve prior to the 1994 LBV event. 



Table 2. Final 'best-fit' values of the model parameters of the 
stars in the HD 5980 binary system for / - 86° and prior to the 
1994 LBV outburst. 





Primary eclipse 




Secondary eclipse 






Primary eclipse 


Secondary eclipse 




Ingress 


Inj 


;ress 


Eg 


ress 




Ingress 


Ingress & Hgress 




.00935 ± .00026 


.00832 


± .00041 


.00481 


± .00024 


star A 


pi = 0.158 


± 0.005 


B4 


.01801 ±.00048 


.01542 


± .00065 


.00922 


± .00043 






Li = 0.398 ±0.021 


Be 


.02604 ± .00067 


.02164 


+ .00080 


.01326 


± .00059 


star B 


= 0.108 


± 0.003 


Bg 


.03351 ± .00083 


.02721 


± .00088 


.01699 


± .00071 




po = 0.269 


±0.014 


Bio 


.04046 ± .00097 


.03226 


± .00094 


.02044 


± .00082 




L2 = 0.300 ± 0.016 





The second step in the procedure consists of solving the 
equations for / - 86°, each half-eclipse being treated in a com- 
pletely independent manner We recall that all radii are reduced 
to the semi-major axis of the relative orbit. 

We first consider the descending branch of the primary 
eclipse. Solutions for pi,p2e,p3e,L2, andjj, are searched for dis- 
crete values of the parameter U2 (0.1, 0.2, 0.3, ... 1) with rather 
broad variation ranges allowed to the variable parameters. Each 
search sequence is based on a series of 1000 trials. A first set of 
results is obtained to provide convergence rates (i.e., an estimate 
of the likelihood of the derived solutions) and a cr value for each 
parameter. A second iteration with reduced variation ranges (the 
original one ±cr) for the parameters leads to a second set of solu- 
tions with improved convergence rates and lower cr values. The 
process is repeated five times, until a stabilization of the param- 
eter values becomes noticeable. For the last iteration, the grid 
of M2 values is enlarged significantly and the respective conver- 
gence rates of the corresponding solutions are used to determine 
the best fit U2 value and to estimate the error on this parameter. 
The obtained results arepi = 0.150 + 0.004,p2e = 0.257 ±0.018, 
P3e = 0.110 ± 0.005, L2 = 0.300 + 0.016, M2 = 0.58 + 0.07, and 
=0.19 + 0.03. 

For the secondary eclipse, solutions are searched for pi , p2a, 
pj,a, L\, and ja- The choice of the parameter u\ for the limb- 
drakening law (Klinglesmith & Sobieski ,1970) of star A has lit- 
tle impact on the other parameters. We repeated the fitting pro- 
cedure for different values of ui (0.1, 0.3, 0.5, 0.7) and recov- 
ered the same solutions within the errorbars. The largest sensi- 
tivity was found forp2n when ui - 0.7. In this rather unlikely 
case, p2a exceeds its usual value by 1 .5 cr. In the following, we 
thus focus on the results obtained with u\ - 0.3, which seems 
a reasonable value for star A (Klinglesmith & Sobieski il970t . 
The descending and ascending branches are analysed separately 
and solutions are searched for discrete values of the parameter 
V (0.1, 0.2, 0.3, ... 1). The same iterative procedure as described 
above for the primary eclipse is applied. The resulting solutions 
for the secondary descending branch are p\ - 0.160 + 0.006, 
P2a = 0.259 + 0.018, p3„ = 0.110 + 0.005, U = 0.385 + 0.028, 
ja = 0.19 + 0.03, and v = 0.60 + 0.20; and for the secondary 
ascending branch pi = 0.163 + 0.008, p2« = 0.290 + 0.019, 
P3fl = 0.105 + 0.003, Li = 0.412 + 0.026, y,, = 0.20 + 0.03, and 
1^ = 0.45+0.15 

The diff'erences between the solutions provided by the three 
half-minima are relatively small compared to the errors and aver- 
age values are derived for the parameters. A distinction between 
absorption and emission values does not appear to be necessary 
any longer. The error in the mean for each parameter p, is com- 
puted in a conservative approach using the expression 



U2 = 0.58 ± 0.07 



v = 0.52 ±0.14 



cr = 



i T:uiPi-<p>f+T:u <T] 

nN 



,2\l/2 



(33) 



where n is the number of values used in the mean and N the num- 
ber of independent sets (half-eclipses) of values. For an analysis 
completed in this way, the values adopted for the parameters of 
the stellar components in the HD 5980 binary system are given 
in Table|2] 



4.3. Discussion 

Figure |4] shows how the synthetic light curve derived from the 
above elements fits the observational data for both the primary 
and secondary eclipses. While the bottom and the wings are 
fairly well fitted, the transparency law that we adopted for the 
extended envelope of star B is probably still too crude to allow 
a perfect match to the observations of the descending and as- 
cending sides of both minima. As a consequence, the size of this 
envelope is probably slightly overestimated by our model. An 
asymmetry in this envelope, inferred by BP91, is difficult to as- 
certain from the present study. The difference between the values 
of p2a provided by the two branches of the secondary eclipse is 
indeed marginally significant only given the errors. L\ +L2 lis 
independent confirmation of an unresolved source of third light 
along the line of sight that accounts for an additional relative 
luminosity of L3 = 0.302. In principle, one would expect a su- 
perior control of the uniqueness of the solution if the bright- 
ness ratios L2/L1 and L3/L1 in the optical could be fixed inde- 
pendently of our light curve analysis. However, in the specific 
case of HD 5980, it is impossible to infer these brightness ra- 
tios from spectroscopy. In this system, there have been changes 
in both the spectroscopic and in photometric properties of the 
binary components, and the brightness ratios are thus epoch de- 
pendent. Spectrophotometric brightness ratios would have to be 
estimated based on the dilution of the emission lines. However, 
as shown by recent spectroscopic studies (e.g. Foellmi et al. 
2008 ), these emission lines contain epoch-dependent contribu- 
tions from components A and B, as well as from the wind-wind 
interaction region. In our analysis of the light curve, we thus 
preferred not to constrain the brightness ratios a priori. Our re- 
sults can however be checked a posteriori against the results 
obtained in the UV domain. Assuming that an 07 supergiant, 
which does not partake in the orbit, was responsible for this third 
light, Koenigsberger et al. (il994j estimated a value of 2.8 for the 
ratio (Li + L2)/L3, which was considered to be consistent with 
the 2.0 obtained from the BP91 elements. It is noticeable that the 
value of 2.3 derived from the present analysis is in even closer 
agreement with the estimate of Koenigsberger et al. (|19941 l. 

Adopting My - -7.5 for the global absolute visual mag- 
nitude of HD5980 (cf. Koenigsberger et al. 1 994 ), the mag- 
nitudes of the binary components are My{starA) = -6.50 and 
My{starE) - -6.19. Because of the difficulties in assigning the 
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Fig. 4. Top panel: the observed light curve of HD 5980 in the Stroemgren v filter as a function of orbital phase. The lower left and 
lower right panels show the synthetic light curves for the primary and secondary eclipses respectively compared to the actual data. 
The synthetic light curves were computed using the mean parameters of the solutions found by our programme. 



various lines of the spectrum of HD 5980 to a specific compo- 
nent, there have been few attempts to establish a full SB2 orbital 
solution (Niemela et al. 119971 Foellmi et al. I2008I I. Therefore, 
our knowledge of the component masses of the eclipsing binary 
remains uncomfortably poor, even though we are dealing with an 
eclipsing system. Foellmi et al. (2008) used a multi-component 
fit of the Niv /1 4058 and Nv /1 4603 lines to derive absolute 
masses of 58 - 79 Mq for star A and 51 - 67 Mq for star B. 
These values differ significantly from the estimates of Niemela 
et al. d 19971 1, who attributed the entire Niv A 4058 line to star A 
and the entire Nv A 4603 emission to star B when inferring mini- 
mum masses of m sin^ / - 28 and 50 Mq, respectively. Adopting 
the solution of Foellmi et al. (I2008I I. the sum of the masses of 
the components of the eclipsing binary equals 109 - 146 M©. 
This corresponds to a semi-major axis of 0.67 - 0.74 AU (144 - 
159 Ro) for the binary orbit, and to stellar radii of 22.7 - 25. 1 R© 



for star A, 15.6 - 17.2 Rq for the core of star B, and 38.7 - 42.8 
R0 for its envelope. 

Alternatively, one could estimate absolute masses from the 
visual magnitudes evaluated above by means of mass-luminosity 
relations derived e.g., from massive star evolutionary models. 
However, this approach is hampered by a lack of precise knowl- 
edge of the spectral types of the stars (and hence their bolometric 
corrections) and by the fact that the components of HD 5980 can 
probably not be considered as 'normal' early-type stars. 

We considered the correlations between the various free pa- 
rameters by plotting them in pairs (see Fig.|5]l. For this purpose, 
we used the results of 10000 trials for each eclipse. In most 
cases, we do not observe obvious correlations; the solutions are 
scattered over a limited part of the parameter plane (see e.g., the 
radius of the core of star B, p3, versus that of star A, pi, de- 
rived from the primary eclipse). However, the solutions show a 
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0.4 



I' I Pi Paa 



Fig. 5. Distribution of the parameters derived from the inversion of the moments of the Hght curve in various parameter planes. 
The left and middle panels are derived from the primary eclipse, whilst the rightmost panel presents the solutions of the secondary 
eclipse. The final solutions derived from each eclipse as well as their 1 - o" error bars are overplotted in each panel. 



clear trend, if we plot the luminosity of star A {L\) as a function 
of the radius of the envelope of star B {p2a, evaluated from the 
secondary eclipse), and fall even along a clearly defined locus if 
we plot the luminosity of star B (L2) as a function of the radius 
of star A (pi, derived from the primary eclipse). These trends 
can be understood at least qualitatively. In fact, for the primary 
eclipse (the opaque core of star A occulting star B), an increase 
in the radius of star A implies that a larger fraction of the light of 
star B is removed. To account for the observed depth of the light 
curve, the total luminosity of star B must decrease. On the other 
hand, during secondary eclipse (star B in front of star A), an in- 
crease in the radius of the semi-transparent envelope of star B 
will produce deeper and broader wings of the secondary eclipse. 
To account for the observed eclipse shape, the model must react 
in terms of an increase in the surface brightness (and hence the 
luminosity) of star A. 

In principle, the photometric variability of HD 5980 as de- 
tected through medium-band filter observations could be af- 
fected by the Doppler shift of strong emission lines that fall 
in the wavelength range covered by the photometric filters. To 
quantify this effect, we simulated an observation of the star 
through the ESO Stroemgren filters. For this purpose, we used 
the spectrum of HD 5980 taken from the spectrophotometric cat- 
alogue of Morris et al. ( 1993) kindly provided to us by Dr J.-M. 
Vreux. Our calculations indicate that less than 1% of the flux 
in the v band comes from emission lines. Therefore, we do not 
expect the Doppler shift to have any significant eff'ect on the pho- 
tometric variability discussed here. The situation would be quite 
different, if we were dealing with data taken in the b band. There, 
about 30% of the flux is produced by emission lines (especially 
the strong Hen /1 4686 line). In this case, the Doppler motion as 
well as variations in the line intensity caused by the wind-wind 
interaction (e.g., Breysacher et al. 119821 Breysacher & Francois 
I2OOOI 1 could lead to significant variations in the observed b mag- 
nitude. 

Finally, we note that the data points of the light curve suggest 
an increase in the brightness of the HD 5980 system at orbital 
phases after (p ~ 0.04 (see Fig.|4]i. Whilst this increase could 
be caused by light reflections from the wind interaction region 
(which should have its concavity roughly turned towards the ob- 
server) near periastron, we caution that this apparent trend is ac- 
tually inferred from observations from a single campaign (cycle 
39 in Fig.|6]l and could therefore be related to intrinsic variabil- 
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Fig. 6. Observed v magnitudes of HD 5980 in the phase interval 
0.04 - 0.16. The various symbols indicate data from different 
orbital cycles. Periastron passage occurs at phase 0.061. 



ity of the WR star rather than represent a genuine phase-locked 
effect. 



5. Conclusions 

We have presented a method that allows us to treat the problem 
of atmospheric eclipses in the light curves of moderately wide, 
eccentric Wolf-Rayet + O binary systems in a semi-analytical 
way. We have then applied this method to the light curve of the 
peculiar system HD 5980 prior to its 1994 outburst. Despite the 
non-uniform sampling of the light curve, the fact that we are 
analysing data from different orbital cycles (hence affected dif- 
ferently by the intrinsic variability of the WR star) and the limi- 
tations of our assumptions on the properties of the WR envelope, 
our method yields consistent results when applied to the primary 
or secondary eclipse. We have been able to constrain some of the 
physical parameters of this system, although the lack of a consis- 
tent SB2 spectroscopic orbital solution prevents us from obtain- 
ing fully model independent parameters. As a next step, we will 
try to generalize our method to the analysis of the eclipses of a 
system harbouring two stars, both with extended atmospheres. 
This should allow us to analyse the light curve of HD 5980 ob- 
served after the LBV eruption. 
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Appendix A: Analytical expressions for the 
moments of the eclipses 

A.1. Primary eclipse (OB star in front) 

Here the relevant radii are: pi, the radius of the OB-star; p2e 
and Pit,, the radii of the W-R core and envelope in emission. We 
note that we have assumed that p2e < pi- Adopting a surface 
brightness law for the W-R star given by Eqs. ( l22l i and ( 1231 ) and 
F{s) = 1 (for s < pi), we find that its total luminosity is given 
by 



/,„ — Lw-R 



(A.l) 



(A.2) 



and 



W-R ■ 



(A3) 



where g(p\,p2e) - min^l,^^) j. The expression of P reduces 



to 

P - Lw-R if pie < pi, 
and finally, 
R„, — 1. 



(A.4) 



(A.5) 



A.2. Secondary eclipse (W-R star in front) 

Here p2a and p^a are the radii of the W-R core and envelope in 
absorption, respectively. Again, we have assumed that p2a < Pi- 
Adopting a transparency law for the W-R star given by Eqs. ( |20] i 
and we find that 



Lo^J(0)np] 
Lo 



1 - Ml + Y 



and 

P = 



1 — Ml 



2miot M^m' 

+ — 

m + I m + 2 ^ 



(A.6) 
(A.7) 



Lo 



\ - U] + ^ 



-«1 -f + 

Pi ^ p\ ) 



+ yah 



\ — U]h + 



Mj h 



.4 3 „6 
1 2u\h 



(A.8) 



where g{p\,p2a) = min^l,(^)^j and h(pi,P2a) 



mm 

Finally, we obtain 



\p2a I V m+lJ 



(A.9) 
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